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Abstract 

The decay rate of a triaxially-deformed proton emitter is calculated in a particle-rotor model, 
which is based on a deformed Woods-Saxon potential and includes a deformed spin-orbit interac- 
tion. The wave function of the / = 7/2^ ground state of the deformed proton emitter ^^^Ho is 
obtained in the adiabatic limit, and a Green's function technique is used to calculate the decay 
rate and branching ratio to the first excited 2+ state of the daughter nucleus. Only for values of 
the triaxial angle 7 < 5° is good agreement obtained for both the total decay rate and the 2"*" 
branching ratio. 

PACS numbers: 21.10.Tg, 23.50. +z, 27.60. +j 
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I. INTRODUCTION 



The fundamental simplicity of the proton decay process in nuclides whose ground states 
are unstable to proton emission has enabled a good deal of nuclear structure information 
to be obtained on nuclei beyond the proton drip line The observable quantities are the 
proton energies and half-lives. In the rare-earth region, the proton emitters are predicted to 
have large static quadrupole deformations 2]. For these cases, analysis of the measurements 



has been carried out using a particle-rotor model, with t 



with an axially-symmetric deformed core 



iyjj,fly,y,Hiiii. The 



le unbound proton interacting 



results of such 



analyses over the past several years has been to obtain a good description of the ground- 
state properties of deformed rare-earth proton emitters, including deformations, occupation 
factors, Nilsson orbitals for the decaying protons, and wave function decompositions. 

Adding to the information provided by the observation of decay protons, recent mea- 
surements have been made of the level structure of the deformed proton emitter ^^^Ho by 
means of in-beam 7-ray spectroscopy Particle- rotor calculations of the energy levels in 
the rotational band lying above the 1 = 1/2" ground state suggest that better agreement 
with experiment would be obtained if the nuclear shape possessed a small amount of static 
triaxial deformation 0|. From the standpoint of proton radioactivity it is therefore of inter- 
est to investigate the effect of a static triaxial deformation on the decay rate of a deformed 
proton emitter. In this work we present such an analysis, and obtain numerical results for 
the decay rate of the deformed proton emitter ^^^^Ho. 

II. COUPLED EQUATIONS IN THE R REPRESENTATION 

We generalize the treatment of Esbensen and Davids [llL extending it to include the case 
of nuclei without axial symmetry. Using Eq. 5A-2 of ref. lls!], we write the wave function of 
an odd-A even-N nucleus consisting of a proton coupled to an even-even triaxially-deformed 
rotor, in the laboratory (space-fixed) system as 



^jM{r,uj) = Y: ^^^^^IIjRtIM), (1) 

IjRr ^ 



where / and j are the orbital and total angular momentum of the particle, R and r are 
the rotational quantum numbers of the rotor, and / is the total angular momentum of the 
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nucleus (I = j + R). The ket \ljRTlM) describes the dependence on spin and angular 
coordinates of the particle and the orientation angle of the rotor, and is given by 



\ljRTlM) = ^ {jmRMR\IM)\RTMR)\ljm). (2) 

mAlfj 

This is the laboratory frame or R- representation as described in Eq. (2) of ref. The 
total Hamiltonian of the proton-core system, 

H = T + V{r, u) + VUr, u) + Hn, (3) 

consists of the relative kinetic energy T, the nuclear plus Coulomb interaction V{t, u), which 
depends on the position r of the proton and the orientation u of the rotor in the space-fixed 
system, the deformed spin-orbit potential Vis{r,uj), and the Hamiltonian if/j of the rotor. 
The detailed parametrization of the nuclear and Coulomb interactions are given in Appendix 
A, and the deformed spin-orbit term is discussed in Appendix B. To proceed, we first expand 
the potential V{r, uj) in the D-functions, which are related to the spherical harmonics 

V{r,u;) = Y.Vx,{r)D',oiO'A'), (4) 

0\ _|_ 1 rl /•27r 

Vx,{r) = j_^d{cose') d<j>'D'^*{e',<P')V{r,e',<P'), (5) 

where 6\ (j)' refer to the angles of the particle with respect to the 3-axis of the rotor. Because 
of reflection symmetry, A and /i are restricted to even values (see Appendix A). Projecting 
with lljRrlM) on the Schrodinger equation H^im = E'^im we obtain a set of coupled 
equations in the radial wave functions 

{h,+Er,r-E)ct>{,R,{r) = - j: H {lJRrIM\Dl,{9\ct>')\l'J'^ir'IM)V^,{r)<^i^,^,Ar) 

I'j'R'r' A>0,M 

- Yl {ljRrIM\VUr,u;)\l'j'R'T'IM)4^,j,,^,{r), (6) 

I'j'R'r' 

where 

and Vo(r) is the monopole part of the Coulomb plus nuclear potential. Here /xq is the proton 
reduced mass, and E^^ is the energy of the rotational state \RtMji). In Appendix B we will 
extract the monopole part of the spin-orbit potential for inclusion in hij. 
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III. THE K REPRESENTATION 



The matrix elements on the ri ght hand side of Eq. © are easiest to evaluate if we go over 



into the K-representation of ref. ll|, which is oriented in the body-centered coordinates of 
the rotor. In this system the quantum number r is identified with the projection Kji of R 
on the rotor's 3-axis. For clarity, the quantum numbers m, Mr, and M in Eq. Q all refer 
to angular momentum projections on the z-axis of the space-fixed coordinate system due to 
the particle, rotor, 3-axis in the body-fixed system are denoted by f2, Kr, and K (see Fig. 
(^). As a consequence, the following relations hold: 

M = Mn + m (7) 
K = Kr + VI. (8) 



In the K-representation we can write the transformed rotor wave function as (see Eq. 4-7 
of ref. 0) 

/ op _|_ 1 

{uj\RTMn) = {u\RKnMR) = ^—-^Df,^j,^{u), (9) 

which is a function of the orientation uj of the rotor in the laboratory frame. For an axially- 
symmetric rotor, we would have Kr = 0, resulting in Q = K from Eq. (jH}. The particle 
wave function is 

\ljm)=Y.DlMm), (10) 
n 

where the single-particle wave function is evaluated in the body-fixed frame of the 

rotor. 

After inserting Eq. Q and ()10|) into Eq. (0), we may contract it using Eq. (lA-43) of 
with the result 



\ljRKnIM) = J^^Y.iJ^RKn\IK)Di,j,iu;)m). (11) 
* °^ K,n 

Because the rotor possesses symmetry after rotating 180° around any of its three axes, 
it is convenient to have the projection K appear only as a positive number. Symmetry 
properties of the wave function require that the quantity Kn = K — Q he an even integer 
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(0, ±2, ±4 . . .) [JJ. We then have: 



(12) 

where Q stands for —Q. Since we are only interested in the low-lying states of the rotor 
where = 0, 2, 4 . . ., we rewrite Eq. ()12|) as 

\ljRKnIM)= Y: A'^^RK^mKIM), (13) 
where 



/2i?+ 1 



A'Ark, = ^i^{3^RKr\IK)^1 + (-1)^ (14) 



and 



\mKiM) = sj^^ [D{,^{u:)m) + {-iy-wi,_^{uj)m)\ . (15) 

For an axially-symmetric nucleus, Kr vanis hes, making Vl = and Eq. ()14|1 and 
become identical to Eq. (14) and (15) of ref. |ll| . 

Inserting Eq. ()13|) into Eq. (1) we can now express the total wave function in terms of 
the new basis flT^ : 



^/M = E E E ^^^mKIM), (16) 
where the radial wave functions are 



ij K>o n 



= E A'n,RKARK,{r). (17) 

RKr 

Note that the triaxial radial wave functions depend on the particle quantum number Q in 
addition to Ij. 

It is easy to show that the amplitudes (fTlj) form an orthonormal transformation between 
the K and the R representation, i.e. 

E 'l2^jn,RKRA^jn,R'Ki,, = ^R,R'^KR,K^n E A^jn,RKR-A^jn' ,RKr = ^n,n'SK,K'- (18) 

K>0 n RKr 

Thus we can transform the results obtained in one representation into the other. After 
inverting Eq. ()17|1 we obtain 



^) = EE^S«i..0SW- (19) 



^IjRKr 

K>o n 
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A. Coupled Equations in the K Representation 

We continue with the evaluation of the matrix elements on the RHS of Eq. Inserting 
expression (jl3p for the spin-angular wave functions in the R-representation we obtain 

{ljRKRlM\D^^,{9', (P')\l'fR'Kn'IM) = 

E E E^Si?/.«(^j'^^^^l^;o(^^0')IO-'f^'^'/M)4^:^,^^, 

K>0 K'>OQ,n' 

and a similar expression for the matrix element of the spin-orbit potential. The advantage 
of using the K representation, Eq. (fT3jl . now becomes evident because each matrix element 
is the product of two parts. The first part involves integrating a product of orthogonal 
D-functions over the orientation coordinates u of the rotor, and yields the important result 
Sk,k'- The second part involves the single-particle wave functions in the body-fixed rest 
frame of the rotor: 

{ijn\D'^,{9',<p')\i'fn') = (-i)^(jiAo|/i)(/r]'A/.|jn). (20) 

Thus 

{ljRKnIM\D'^,{d',<p')\l'fR'KnJM) = T.^'n'RKjmDU(^'AWmAfn',R'K,r 

K'>OQ,n' 

(21) 

For the spin-orbit potential we have a similar expression. We show how to evaluate the 
matrix elements of the spin-orbit interaction in Appendix B. Because the matrix elements 
are diagonal in K, Eq. (jSI) shows that the particle projections Q, Q' are restricted by 1^7 — ^2' | 
= an even integer. We still have the previously determined restrictions A even and |ir — il| 
= an even integer. 

We now obtain the coupled equations in the K representation by multiplying Eq. © by 
^jn,RKii ^iid summing over RKr, using Eq. (fTTjl and (fT8|l . The presence of the rotational 
energy Ejikj^ requires the use of Eq. (fT^ for the radial wave function. Thus we obtain 



K'>o n' 



EE E m\DU(^',<i>'wmvx,ir) + mivMii'm \ ^\%{t\ (22) 
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where 

WjUQ' = H ^jn,RKiiERKaAjn> ^rKr- (23) 

RKr 

To recover the axially-symmetric case we set Q = Q' = K, and /i = in Eq. ()22j) and ()23p. 
This removes the sums over fi', i^^^, and /i. 

B. Adiabatic limit 

rn 

For a comparison with the results obtained in 11], we will solve the coupled equations 
fl22|) in the adiabatic limit, where the rotational energies Erkr of the core are set to zero. 
This sets equal to zero the second term on the left-hand side of Eq. (I22j) . In this case we 
have 

I'j'W [\>o,^l ) 

(24) 

The important thing to notice is that the coupled equations are then diagonal in K, but 
include fi- mixing. In other words, in the adiabatic limit, K is still a good quantum number, 
but, in addition to the Q = K component, the interaction mixes into the wave function 
components with Q = K ± 2, K ± 4, etc, subject to the restriction \Q\ < j. It should also 
be noted that the Clebsch-Gordan coefficient in Eq. ()20|1 is non-zero only when Q' + fi = Q. 

The consequence of this f2-mixing is that the number of coupled equations will be sub- 
stantially larger than is found in the axially-symmetric case. As an example, consider the 
combinations of j and Q needed to solve the equations for the I = K = ^ ground state of 
the deformed proton emitter ^"^^Ho. The spherical states involved will be j = | , | , ^ ; 
Y , and y . In addition to fl = K, for each j there will be associated up to 7 more values 
of Q. Table d shows the permissible j, fl combinations for the I = K = ^ ground state 
of ^^^Ho. The total number of wave function combinations is 30, which is to be contrasted 
with only 5 for the axially-symmetric case. 

C. Decay Rate Calculation 

We obtain the partial decay rate for proton emission from a state having angular mo- 
mentum J = in the adiabatic limit via either the direct method (Dir) or the distorted 



7 



wave Green's function method (DW) using Eq. (7) of 11|. While in the axially-symmetric 
case the daughter states were labeled only by the quantum number R, in the triaxial case 
we need the second label r or Kji as well: 

^ RKr - J- IjRKn - —— 2^ V^ljRKu I ) l^^J 



l3 ^ Ij 



where 



^mZ = -Wr ^ r'drFi{knr){ljRK^IM\V{r,u;) + VUr,uj) 

-^\l'j'R'Kn'IM)4^,^,^Jr). (26) 

Here is the proton wave number for the decay to the daughter state with quantum 
numbers RKr- Having obtained the radial wave function (pij^ir) for a given K = I in the 
adiabatic limit, we can construct the associated radial wave functions in the R representation 

'/'U«(0 = E4S^/..</>SW (27) 
n 

by means of Eq. (|T9|) . Inserting this wave function into Eq. (j26|) . we obtain for the Direct 
method 

Using this equation we can immediately see the effect of the triaxial shape on the decay 
rate: for decay to the ground state {R = 0), the Clebsch-Gordan coefficient has the value 
Sj,i^n,K- Thus only the Q = K component of the wave function will participate in this decay 
branch, effectively reducing its decay rate relative to the that for the axially-symmetric case 
since there are now other fi-components present in the triaxial wave function. 

For the Green's function method, after performing the same operations that led to Eq. 
fl21|) to obtain the matrix elements of the interaction in Eq. ()26|) . the second orthogonality 
relation in Eq. (fTH)) allows the summations over R' and Kjii to be performed. The final 
expression for the Green's function method is 

4S = -^T.T.Ajn,RK, r'drF,{knr){mV{r,u;) + VUr,uj) 
'^R I'j' n,Q' •'^ 

-^\l'm4rn'ir). (29) 



This expression reduces to the analogous one for an axially-symmetric nucleus when K is 
set equal to Vt and the sums over /i, VL and VL' are removed. 



IV. APPLICATIONS 

In this section we use the formalism described above to calculate decay rates for the 
deformed proton emitter ^^^Ho (/ = = 7/2~), in both the axially-symmetric and triaxial 
cases. The calculations have been performed in the adiabatic limit, using the potential 



parameters found in ref. The ground-state wave function was expanded in spherical 

components with j = 7/2~, 9/2^, 11/2", 13/2", and 15/2". All multipole expansions were 
carried out up to A = 14. The results are compared with experimental values for the total 
decay width and the branching ratio for decay to the first excited 2"*^ state in the daughter 
nucleus ^^°Dy. 

A. Axially-Symmetric Case: the Deformed Proton Emitter ^^^Ho 

To verify the formulation, we have computed the decay half-life and branching ratio to 
the 2"*" state in the daughter nucleus for the deformed proton emitter ^^^Ho , using deformed 
spin-orbit matrix elements for the axially-symmetric case (7 = /i = 0, Q = Q' = K) 
calculated from Eq. ()B17|) . which, again for A > 0, now reads 

{ljm,{rMl'fK) = 2KoE ^^^^S=^(^JII [V/A,(r)Dio(f)]. (-.V x (30) 

A>0 V + i L 

where the reduced matrix element is that calculated in Eq. ()B16|) . The results of the 
calculations are shown in Table IHl and are identical with those obtained using the axially- 
symmetric deformed spin-orbit term of ref. |ll|. It should be noted that the numbers 



obtained here and in ref. 



ll| differ slightly, since the 2"*" excitation ener gy of 202 keV in 



known when the calculations in ref. 
of 166 keV was used. 



the daughter nuclide ^"^^Dy used here has only recently been measured 0, Q8| , and was not 



111 ] were carried out. In that work an excitation energy 
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B. Triaxial Case: the Deformed Proton Emitter Ho 

As previously mentioned, particle-rotor calculations of the energy levels in the rotational 
band lying above the / = 7/2~ ground state of the deformed proton emitter ^^^Ho were 
reported in ref. These calculations suggest that better agreement with experiment 

would be obtained with the introduction of a small amount of static triaxial deformation. 
In order to assess the effect on the proton decay rates of adding such a deformation, we have 
used the formalism developed in the present work to compute, in the adiabatic limit, the 
decay rate and 2^ branching ratio for this proton emitting nuclide, as a function of the angle 
of triaxiality 7. The calculations were performed with (32 = 0.244, /?4 = —0.046. Preliminary 
results have been reported in ref. jiol . 

Fig. 121 shows the product of calculated total decay width Tcaic and spectroscopic factor 
Scale = V? = 0.6 for the ^^^Ho ground state, plotted as a function of the triaxial angle 7. 
The spectroscopic factor was obtained from a BCS calculation, using a proton pairing gap 
Ap of 0.9 MeV. The shaded area represents the experimental measurement, and the error 
bar attached to the calculated curve represents the uncertainty in the calculated width due 
to the uncertainty in the proton energy. It is seen that calculation and experiment agree 
well for small values of 7. 

Fig. El shows the calculated 2"*" branching ratio, plotted as a function of 7, along with the 
experimental value, 0.0070(15) The small error bar on the calculated curve is due to 
the uncertainty in the proton energy. The agreement between calculation and experiment is 
excellent for small values of 7, and this suggests that a static triaxial deformation, if present, 
is limited to an angle of < 5°. 

Plotted in Fig. El are the amplitudes for the allowed fi-values of +7/2^, +3/2^, —1/2^, 
and —5/2^ for the j = 7/2 spherical component of the ^^^Ho ground state wave function. 
As expected, for 7 = 0, only the Q = +7/2~ component is present, but with increasing 
7, other ^-components differing by ±2, ±4 . . . begin to mix into the wave function. The 
decay proceeds primarily to the daughter ground state, and only the Q = +7/2^ component 
participates in this branch. Thus the decrease in the total decay width with increasing 7 
seen in Fig. |21 tracks with the decrease in the Q = +7/2^ wave function component, which 
itself follows from the increasing appearance of other fi-values in the wave function. 
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V. CONCLUSIONS 



In this work we have developed a formahsm to include the effect of static triaxial defor- 
mation on calculations of the decay rate for a deformed nucleus. The main complications 
over the axially-symmetric case are the additional dependence of the interaction potential 
on the proton azimuthal angle 0, and a consequent increase in the complexity of the de- 
formed spin-orbit potential. The extra dimension causes additional (^-components to be 
introduced into the wave function. The matrix elements of the spin-orbit interaction have 
been calculated using a tensor algebra approach, as described in Appendix B. 

We have applied the abovementioned methodology to calculations of the decay rate of the 
deformed proton emitter ^^^^Ho, which has spin-parity 7/2^. In ref. [ll] it was shown that 
the non-adiabatic coupled- channels approach does not yield results for ^'^^^Ho decay which 
agree with experiment, either for the absolute decay rate or the branching ratio for decay to 
the first 2"^ of the daughter nucleus ^^°Dy. This is because of the presence of Coriolis mixing, 
which is particularly strong for high spin states. Empirically it is known that it is necessary 



to quench the Coriolis mixing in order to obtain good agreement with data [20|, |21[. It 
was found in that good agreement with experiment was obtained in the adiabatic limit, 
where the energies of the rotational states of the daughter nucleus are set to zero. For this 
reason we have performed the triaxial calculations in the adiabatic limit. 

After first checking the results of the calculation against an axially-symmetric code, tri- 
axiality was introduced, ranging up to a 7 angle of 40°. While the sensitivity of the resulting 
total decay rate to triaxial angle 7 was not high, the calculation of the branching ratio for 
decay to the first 2+ state of the daughter showed a strong dependence on 7, essentially 
ruling out angles greater than 5°. We believe that the branching ratio calculation is quite 
reliable, since factors such as absolute spectroscopic factors tend to cancel. 

In conclusion, we do not believe that a static triaxial deformation plays an important 
role in helping to explain the decay rate of the deformed proton emitter ^^^^Ho. Triaxiality 
may still exert an influence on the nuclear structure of this nuclide, but most likely at higher 
spins and excitation energies. 
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APPENDIX A: PARAMETRIZATION OF THE INTERACTION POTENTIAL 



For triaxia^ 
Hill-Wheeler 



nuclei, the degree of triaxiality is denoted by the angle 7, which is one of the 
22| | coordinates. We parametrize the nuclear interaction between the valence 



proton and the deformed core nucleus in terms of the Fermi function f{x) = [1 + exp{x)] ^ 



as 



V'„(r.».^) = </(^^^). (Al) 



where vj^^ is the depth of the nuclear potential, a{6) is an angle-dependent diffuseness as 



in Eq. (A6) of and 

{2 4 
1+ T.^2Md,<P)+ T.^'^M^^^)}- (A2) 
^=—2,even ^=—i,even 

Here {9, 0) are the angles between r and the 3-axis of the core, and 

(32 sin7 



020 = 1^2 COS7 022 = (32-2 



^40 = |/^4 (5 COS^7 + 1) 

a42 = a4_2 = — 31^/94 "\/30sin27 

044 = a4_4 = ^/?4 v^TO sin^7. (A3) 
The radius is calculated as ^ 

RN = rJ , 1 ' , (A4) 



,C(/32,/54,7), 

where 

C(/32, /34) = / P I 1 + E a2,Y2^{e, 0) + «4^>'4m(^, 0) I (A5) 



/ ^ 1 + E <^2,Y2,{d, 0) + «4^^4; 
\ ^=—2,even ^=—4,even 



is the volume preserving factor and is the mass number of the core. Equations ()A2|) 
and ()A5|1 reduce to Eq. (Al) and (A2) of ^11^ for the axially-symmetric case (/x = 7 = 0). 
Values of tq = 1.25 fm and a = 0.65 fm are used in our calculations. 

We para,„etn.e the charge deus.ty of the eore a similar way. usurg thes.ight.y different 
radial parameters of Tq = 1.22 fm and ac = 0.56 fm, as was done in ref. |11| . 
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APPENDIX B: DEFORMED SPIN-ORBIT POTENTIAL 



For the solution of the coupled equations Eq. (j24j) . matrix elements of the spin-orbit 
potential {ljQ\Vis{r,6,(f))\l'j'Q') are needed. As in ref. we use a deformed spin-orbit 
potential, but with the addition of (/)-dependent terms not found there. We follow the 



tensor algebra approach of Hagino 



23l |. who has calculated matrix elements of the spin- 



orbit interaction in the context of particle-vibration coupling for spherical nuclei. 

The starting point is to express the deformed spin-orbit interaction in the Thomas form 
111 ] with the deformed Fermi function of Eq. ()A1|) : 



ViJr,9, 



4K, 



V/ 



xp-s 



We now make a multipole expansion of the Fermi function 



2A -1-1 r 

/amW = I Di;{r)f{r,9,<p)d{cos9)d<p, 



(Bl) 



(B2) 



dr 



2A 



47r 



^^df{r,e,<P) 



dr 



d{cos9)d(j). 



(B3) 



So 



4Ko(E[V/A,(r)Z^io(f 
'^dfooir 



xp-s 



4K 



' \>0,fj. 



xp-s 



4v;, 



}^dfoo{r] 
r dr 



Is 



I E [V/. 



AAO^io(f) 



X>0,fj. 



'-iV X a) 



(B4) 



The first term on the RHS of Eq. ()B4|) is the monopole part of the spin-orbit potential, 



4K,o- 



^dfQo{r] 



Is, 



(B5) 



r dr 

which can easily be incorporated into the monopole part of the Hamiltonian. The second 
term, 

5Vis{r,9,<P) = Vis{r,9,<P) - Vi'^.ir), (B6) 

can be decomposed into a sum of angular momentum tensors involving the spherical har- 
monics Y\fj_{r:), using the gradient formula Eq. (5.9.17) of {2^: 



V/v(r)Dio(f 



-iV X a) 



y/MA + l) f dh,{r) A , 
2A + 1 \ dr r^'^^\ 



[rA+i(-^V X a)] 



Xfi 
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2A + 1 I dr 



A + 1 



fx,{r) [1a-i(-^V X a)] 



A/1 



1^ {^Qx,{r)Th^ - V^Px,{r)Tl^] , (B7) 



2A + 



where 



Ti^ = [Fa±i(-^V X a)] 



(B8) 



Qx,{r) 



dfxiijr) _ A 

dr 
dfx,{r) 



fx^{r), and 

r 

+ fx^,{r). 



dr r 

To solve the coupled radial equations we need to obtain the matrix elements of the oper- 



TiA/i 



ator J2x>o.fi fxii{r)D^Q{r) ■ (— iV x a), or specifically, the matrix elements of 
between spin-angular momentum states. This is done by using the Wigner-Eckart theorem 



VW 



(B9) 



along with reduced matrix elements given in Eq. (58) and (59) of ref . |25l and Eq. (A2.14) 
of ref. j^. Eq. ()B9|) is the point of departure between axially-symmetric and triaxial 
calculations of the spin-orbit matrix elements, since the reduced matrix elements are inde- 
pendent of the projection quantum numbers Q', Q, and /x. For axial symmetry we will set 
n = n' = K and = n = 0. 
Putting L± = A ± 1 we have 

{lj\\Ti\\l'f) = [Y^.i-tVxafwi'f) 

= 5](-l)^+^+^'+y36(2A + l)(2j' + l)(2j„ + 1) 




/' 1 

1 i 1 

2 2 
Ja f 1 



}mYLjLja){yij'^\\yi'), (bio) 



where 



/ (2j + l)(2L^ + l)(2j„ + l) J L± J, 



47r 



(BU) 



and 



(F,J|V||Fr) = V2F+1(/'010|L0) 



dr r 2r 



(B12) 
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From the properties of the Clebsch-Gordan coefficient in Eq. ()B12|) it can be seen that la 
can take on only the 2 values /' ± 1, indicating that the parity of laja is opposite to that of 
I'j' and Ij. The permissible values of ja in the summation in Eq. ()B10|) are ja = ± 1/2, 
and may be further restricted by the triangle relations for the 6j-symbol, namely A{L±jja) 
and A{j'lja). Additionally, for the calculation of the matrix elements, the d/dr term in 
Eq. ()B12|) needs to be modified, since the coupled equations are in the wave function 0/(r), 
while the Hamiltonian Q acts on the wave function (^, which contains (f)i{r)/r. Since 

( <i)i{r) \ _ d(f)i{r) 4>i{r) 
dr \ r J dr r ' 

we see that the operator d/dr must be replaced by {d/dr — 1/r) in the coupled equations 
for (f)i{r). This replacement has been performed below. 

Finally, after noting that A is always even in this application, we have 



2A + 
where 



lot 3 a. 



^ + ^{ni' +i)-L{L + i)} 

dr 2r 



(B14) 



36(2/ + 1)(2L± + l)(2j + 1)(2/^ + 1) 
(2A + 1) 



(B15) 



and 



AHjaj'jL^) = {2ja + 1) 



L+ 1 A 



3 3 3c 



la V 1 

1 i 1 

2 2 

3' 1 



3 L± Ja 

1 n 1 



(/'010|/^0). 



After some rearrangement, we can write the reduced matrix element of the operator 
V/A,(r)Z^;io(f)] ■ (-^V X a) as 

{lj\\ [v/A,(r)Dio(f)]-HVxa||/y) = 

I VXQx,{r)C\j'jL.) Y: A\ja3'3L-) - VxTlPx,{r)C\j'jL^) ^ A\ja3'3L+) \ 

+ ^ W\Q,,ir)C\fjL.) Y: A\ja3'3L.) [/'(/' + 1) - + 1)] } 

~Yr ¥^P^^^^)C\3 3^^) E A\ja3'3L+) [I'i'l' + 1) - UL + ■ (B16) 



d 
dr 
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Recall that our original goal is to compute the spin-angular coupling matrix elements of 
the deformed spin-orbit interaction. For A > we have, from Eqs. ()B4|) and ()B9|) . 



A>0,M V + ^ 



m Vfx,{r)D'{v).{-tVxa\\l'j'). (B17) 
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TABLE I: Combinations of j and Q for I = K = originating from i = | ~^ ^ ■ "^^^ bullet 



(•) denotes states occurring in the axially-symmetric case. 



j 




Axially-Symmetric 


J 


o 

1 i 


Axially-Symmetric 


7- 
2 


5 
2 




13- 
2 


13 
2 




7- 
2 


1 
2 




13- 
2 


9 

2 




7- 
2 


3 
2 




13- 
2 


5 
2 




7- 
2 


7 
2 


• 


13- 
2 


1 

2 




9- 
2 


9 
2 




13- 
2 


3 
2 




9- 
2 


5 
2 




13- 
2 


7 
2 


• 


9- 
2 


1 
2 




13- 
2 


11 

2 




9- 
2 


3 
2 




15- 
2 


13 
2 




9- 
2 


7 
2 


• 


15- 
2 


9 
2 




11- 
2 


9 
2 




15- 
2 


5 
2 




11- 
2 


5 
2 




15- 
2 


1 
2 




11- 
2 


1 
2 




15- 
2 


3 
2 




11- 
2 


3 
2 




15- 
2 


7 

2 


• 


11- 
2 


7 

2 


• 


15- 
2 


11 
2 




11- 
2 


11 

2 




15- 
2 


15 
2 





TABLE IL Comparison of calculated and experimental values for ^^^110(7/2 ) total decay rate and 
2+ branching ratio. Input values are: £"^(2+) in i^°Dy = 202 keV, P2 = 0.244, and (3^ = -0.046. 
Calculations were done in the adiabatic limit , with no static triaxiality, and uncertainties coming 
from the measured proton energy of 1169(8) keV ,15i] and the decay half-life of 4.2(4) ms have 
been taken into account. The experimental spectroscopic factor S^xp 

is defined as Texp/^caic- The 
calculated spectroscopic factor Scale = is obtained from a BCS calculation (see text). 



r,aic (10-19 MeV) 


Tcxp (10-19 MeV) 


Scale 


Sexp 


BR{2+)calc 


BR{2+)cxp 


i.5il:i 


1.09(10)" 


0.6 




0.0071(5) 


0.0070(15)^ 



''Ref. S 
^Ref. [16 
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FIG. 1: (a) Relationship of the angular momentum vectors j, R, and I in the axially-symmetric 
case, (b) Same as (a) except for the triaxial case. In both cases K is the projection of the total 
angular momentum I on the 3-axis of the body-centered system. 

FIG. 2: Galculated total decay width, in units of 10~^^ MeV, for ^'^^S'Ho, plotted as a function 
of the triaxial angle 7. The shaded area represents the experimental measurement, and the error 
bar attached to the calculated curve represents the uncertainty in the calculated width due to the 
uncertainty in the proton energy. 

FIG. 4: Amplitudes for various O-values of the j = 7/2~ component of the ^^^^Ho wave function, 
plotted as a function of the triaxial angle 7. 



FIG. 3: Calculated branching ratio in percent for the decay of ^^^^Ho to the 2+ state of ^^'^Dy 
at an excitation energy of 202 keV, plotted as a function of the triaxial angle 7. The shaded 
area represents the experimental measurement, and the error bar attached to the calculated curve 
represents the uncertainty in the calculated branching ratio due to the uncertainty in the proton 
energy. 
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^"•^SHo Proton Decay Width vs. Triaxiality 
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